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OBJECf IVE — Growth hormone (GH) deficiency is associated with insulin resistance and 
diabetes. The aim of the current study was to determine incidence of diabetes during GH re- 
placement therapy (GHRT) and the effect of GHRT on fasting plasma glucose concentrations and 
HbA lc in adult patients with GH deficiency. 

RESEARCH DESIGN AND METHODS— A total of 5,143 GH-deficient patients (male 
49.9%; mean age ± SD, 49 ± 13 years; BM1 29.1 ±5.9 kg/m 2 ) were analyzed. Mean observation 
period was 3.9 years (range 0.01-13). Total number of patient-years was 20,106. Observed 
number of cases (O) was compared with expected number of cases (E). Reference rates were 
from Sweden, three additional European regions, and one U.S. region. 

RESULTS — Patients who developed diabetes (n = 523) were older; had higher BMI, waist 
circumference, triglyceride concentrations, and blood pressure; and had lower HDL-cholesterol 
concentrations (P < 0.0001) than those who did not develop diabetes. Diabetes incidence was 
2.6 per 100 patient-years, equal in both sexes, and significantly increased compared with the 
Swedish reference (O/E = 6.02; P < 0.0001) as well as with the four other populations (O/E = 
2. 1 1-5.22). O/E increased with BMI and decreased with duration of GHRT (P < 0.0001). There 
was no significant association with GH dose (P = 0.74) or IGF-I SDS (P = 0.47). In subjects not 
developing diabetes, plasma glucose concentrations increased from 84.4 ±0.9 mg/dL to 89.5 ± 
0.8 mg/dL (0.70 mg/dL/year) and HbA lc increased from 4.74 ± 0.04% to 5.09 ± 0.13% 
(0.036%/year) after 6 years of GHRT. 

CONCLUSIONS — Diabetes incidence appears to be increased in GH-deficient patients re- 
ceiving GHRT and exhibiting an adverse risk profile at baseline. Therefore, glucose homeostasis 
parameters should be monitored carefully in these patients. 
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Growth hormone (GH) plays an im- 
portant and complex role in glu- 
cose, lipid, and protein metabolism 
in humans (1). It decreases glucose oxi- 
dation and glucose uptake by muscle and 
increases gluconeogenesis. In addition, 
GH stimulates lipolysis, lipid oxidation, 
and protein synthesis and decreases pro- 
tein breakdown, amino acid degradation/ 
oxidation, and hepatic urea formation (1). 
GH effects are modified by fasting, food 
intake, stress, and drugs. 

Both GH deficiency and GH excess 
are associated with disturbances of car- 
bohydrate metabolism (2-6). Therefore 
special attention has been paid to changes 
in glucose and HbA lc concentrations as 
well as glucose metabolism in GH-deficient 
patients during GH replacement therapy 
(GHRT). It has been speculated that pos- 
sible adverse GH effects, leading to insulin 
resistance and hyperglycemia, might be 
counterbalanced by the favorable effects 
on body composition, resulting in muscle 
mass increase and fat mass decrease. Earlier 
studies report equivocal effects of GHRT 
on insulin sensitivity in GH-deficient pa- 
tients (7-9). However, in these studies, 
only a small numbers of patients were 
included and the weight-based GH dose 
was higher than the currently used dose 
aimed at normalizing IGF-I. Indeed, it has 
been demonstrated that low-dose GHRT 
over a period of 12 months enhances in- 
sulin sensitivity, whereas GH in a so-called 
standard dose (0.48 mg/day) had no effect 
on glucose metabolism (10). 

Here we report on the incidence of 
diabetes mellitus in a cohort of 5,143 
GH-deficient patients receiving GHRT 
for up to 13 years. In addition, results of 
a longitudinal study looking into yearly 
fasting plasma glucose and HbA lc con- 
centrations in a subgroup of patients over 
an observation period of 6 years are re- 
ported. 
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RESEARCH DESIGN AND 
METHODS— Out of the 13,294 pa- 
tients in KIMS (Pfizer International Met- 
abolic Database, data lock July 2008), 
only patients with severe adult-onset GH 
deficiency confirmed with an accepted 
GH stimulatory test (11) and naive to GH 
treatment were selected. Moreover, pa- 
tients with a history of Cushing disease 
or acromegaly, with diabetes at baseline, 
or who were <20 or >80 years were ex- 
cluded. In addition, in patients with 
yearly control examinations and who 
did not develop diabetes, fasting plasma 
glucose concentrations and HbA lc levels 
were analyzed in a longitudinal study 
spanning an observation period of 6 
years. 

Diabetes was diagnosed according to 
the criteria of the American Diabetes 
Association (12) and/or when reported 
as an adverse event and/or when antidia- 
betes drugs were prescribed. 

Plasma glucose and HbA ic values 
were determined in local laboratories; 
lipid and IGF-I concentrations were de- 
termined in a central laboratory. Serum 
IGF-I measurements were performed at 
Kabi Pharmacia, Stockholm, Sweden, be- 
tween January 1994 and October 1997 
and thereafter at Sahlgrenska University 
Hospital, Gothenburg, Sweden, using a 
radioimmunoassay acid/ethanol precipi- 
tation of IGF-binding proteins (Nichols 
Institute Diagnostic, San Juan Capistrano, 
CA) until November 2002, a chemilumi- 
nescence immunoassay (Nichols Advan- 
tage System) until September 2006, and 
currently, the Immulite 2500 (DPC Sie- 
mens, Eschborn, Germany) (13). IGF-I 
data are presented as SDS. SDS refers to 
the general population and corresponds to 
Z score = (observed value - mean value of 
the reference population)/SD value of 
reference population. The reference 
range ± 2 SDS was calculated from an 
age-corrected algorithm in which IGF-I 
was expressed in micrograms per liter. Be- 
cause the methods and references have 
changed over the years, different SDS for- 
mulas have been used. The SDS algorithm 
formula used between 1994 and 2002 was 

ln(IGF-I) -/in(age) 
SD ln 

The mean function / ln (age) on the natural 
logarithmic scale and the SD was (5.95 - 
0.0197 X age) and 0.282, respectively, 
between 1994 and 1997. Between 1997 
and 2002, the mean function was (5.92 - 
0.0146 X age) and the SD was 0.272. 



From 2002 on, the algorithm of Brabant 
et al. (14) has been used. All central 
measurements were performed on the 
blood samples, which were routinely 
obtained. 

Statistical analyses 

Descriptive statistics are presented as 
mean ± SD or range (depending on the 
skewness of the respective distribution). 
Crude t tests were used for testing differ- 
ence between patients with and without 
diabetes mellitus, respectively. The asso- 
ciation between diabetes occurrence and 
sex, BMI (kg/m 2 ), family history of diabe- 
tes (yes/no), age (years) at diagnosis of 
pituitary disease, age at diagnosis of GH 
deficiency, age at GHRT start, years be- 
tween diagnosis of pituitary disease and 
GHRT start, years since GH treatment 
start, attained age during follow-up, 
number and nature of additional pitui- 
tary deficiencies, surgery before KIMS en- 
try, and calendar period were assessed as 
described below. In a special 1-year left- 
truncated analysis, the association between 
diabetes and IGF-I SDS and treatment GH 
dose at 1 year from treatment start was as- 
sessed. Observed number of cases-to- 
expected number of cases (O/E) ratios of 
patients diagnosed with diabetes were cal- 
culated, stratified for age and sex. The ex- 
pected number of cases was calculated 
using normal population-based age- and 
sex-specific reference incidence rates 
from a recently published study on diabe- 
tes incidence in the Kronoberg county of 
southern Sweden (15) and correspond- 
ingly stratified patient -years in KIMS. In 
addition, O/E ratios were calculated using 
four other reference populations from 
central and southern Europe as well as 
the U.S. (16-19). Patient-years were cal- 
culated from GH treatment start within 
KIMS to date of diabetes, latest known 
visit date, date of death, or KIMS exit 
date, whichever came first. The observed 
number of diabetes cases was assumed to 
follow a Poisson distribution. The 95% 
CIs were calculated using Byar approxi- 
mation formula (20). More detailed mod- 
eling of O/E ratios, incorporating variables 
as above, was performed using Poisson re- 
gression (PROC GENMOD, SAS software 
version 8; SAS, Marlow, U.K.) fitted on 
grouped data. The analyses used the ex- 
pected number of cases as an offset vari- 
able. In a special analysis, glucocorticoid 
replacement in all patients and testoster- 
one replacement in males were analyzed as 
time-dependent covariates. For trend tests, 
category (group) mean values were used 



as numerical representations. Signifi- 
cance level was set at P < 0.05. 

With control for age and sex, estimates 
of the yearly change in plasma glucose 
concentration and change in HbA lc were 
calculated in subgroups of patients (n = 
293 for plasma glucose and n = 451 for 
HbA lc ) with complete information at GH 
treatment start and at the following 6 
yearly visits. Repeated-measurement re- 
gression was applied. Variance-covariance 
matrix was estimated as unstructured. 

RESULTS— Data from 5,143 patients 
(male-to-female ratio 49.9:50. 1 ; mean age 
at KIMS entry ± SD, 49 ± 13 years; BMI 
29.1 ± 5.9 kg/m 2 ) were analyzed. Mean 
follow-up time was 3.9 years (median 3. 1 , 
range 0.01-13), summing up to a total of 
20,106 patient-years. In 649 patients, 
treatment duration was >8 years. Seven 
countries contributed 85% of the patients 
and 81% of the patient-years (Belgium, 
Germany, the Netherlands, Spain, Sweden, 
the U.K., and the U.S.). 

A total of 523 (263 males and 260 
females) of the 5,143 patients (10.2%) 
developed diabetes after a median period 
of 1.7 years (range 0.02-10.3). Baseline 
characteristics of GH-deficient patients 
developing diabetes during GHRT and 
patients not developing diabetes are 
shown in Table 1 . Patients who developed 
diabetes were significantly older; had 
higher BMI, waist circumference, waist-to- 
hip ratio, and triglyceride concentrations; 
and had significantly higher systolic and 
diastolic blood pressure and significantly 
lower HDL-cholesterol concentrations. 
Total and LDL-cholesterol concentrations 
were comparable in the two groups, as 
were GH dose and IGF-I SDS at 1 year of 
GHRT. 

Overall diabetes incidence was 2.6 
per 100 patient-years, gradually decreas- 
ing from the highest value of 4.1 per 100 
patient-years during the 1st year of GHRT 
to 1 .0 per 100 patient -years after > 8 years 
of GHRT (Fig. 1). Incidence was equal in 
both sexes. The overall O/E ratio, com- 
pared with the incidence rate in Kronoberg 
county, was 6.0 (95% CI 5.5-6.6; P < 
0.0001) and decreased from 10.8 in the 
1st year of GHRT to 1.9 after 8 years 
(P trend < 0.0001) (Fig. 1). When the 
diabetes incidence in KIMS was com- 
pared with incidence rates in age- 
matched populations of other European 
regions and the U.S. (16-19), the O/E 
ratios ranged from 2.11 (1.86-2.39) to 
5.22 (4.76-5.72); for more details, see 
Supplementary Table 1 . 
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Table 1 — Characteristics at GHRT start and GH dose and IGF-I SDS at 1 year of 
subjects developing diabetes or not 





Nondiabetic subjects 


Diabetic subjects 


P value 


n 


4,620 


523 




Sex (n) 


Male 


2,301 


263 








260 




Age at KIMS entry (years) 


48.5 (12.6) 


52.7(12.0) 


<0.0001 




Waist circumference (cm) 


Female 


92.7 (14.3) 


99.4 (14.4) 


<0.0001 


Waist-to-hip ratio 


Male 


0.97 (0.08) 


0.99 (0.08) 


<0.0001 


Female 


0.87 (0.09) 


0.89 (0.08) 


<0.0001 


Total cholesterol (mg/dL)* 


227.4 (47.2) 


231.2 (51.0) 


0.19 










HDL (mg/dL) * 


Male 


44.5 (13.9) 


41.0(10.8) 


<0.0001 


Female 


54.9 (15.9) 


48.3 (13.5) 


<0.0001 










Systolic blood pressure (mmHg) 


127.9 (18.0) 


132.9(19.6) 


<0.0001 


Diastolic blood pressure (mmHg) 


79.6 (10.7) 


81.1 (11.4) 


0.005 


GH dose at 1 year (mg/dL) 


0.36 (0.22) 


0.36 (0.24) 


0.83 


IGF-I SDS at 1 year 


0.39 (1.44) 


0.50 (1.43) 


0.21 



Data are mean (SD), except for n. 'Conversion factor 38.67 mmol/L to mg/dL. §Conversion factor 88.57 
mmol/L to mg/dL. 



In the final Poisson regression analysis, 
sex, BMI, attained age during follow-up, 
years between pituitary diagnosis and 
GH treatment start, and years since first 
treatment were statistically significantly 



associated with occurrence of diabetes, 
whereas the other studied variables were 
not (Tables 2 and 3). O/E ratio compared 
with the Swedish population increased 
1.36% per year for year from diagnosis 




I 



I 



I ■ 



I 



25 



20 



15 



O 
3 

10 §■ 



< 1 1-2 2-3 3-4 4-5 5-6 6-7 7-8 >8 

Years 

Figure 1 — Diabetes mellitus by years since the start of GHRT. Diabetes incidence per 100 pa- 
tient-years (gray bar) and O/E cases when compared with a reference population in southern 
Sweden (15). O/E ratio (black bar) P trend < 0.0001. O/E ratio trend estimate in Poisson re- 
gression model -153%/year (95% CI -11.5 to -19.0; P < 0.0001). 



of pituitary disease to KIMS entry 
(95% CI 0-2.66) and decreased 15.3% 
per year for year since treatment start 
( 1 1 . 5 - 1 9 . 0) . For BMI , the O/E ratio in- 
creased 9.4% per unit BMI (7.8-10.95). 
Females had 22.9% higher O/E ratio than 
males (2.6-47.2). O/E decreased by ~4% 
per year of attained age (P < 0.0001) 
with a small quadratic adjustment (Ta- 
ble 3) . Table 2 shows the explanatory value 
for each included variable in the multiple 
regression model. 

No statistically significant association 
between IGF-I SDS at 1 year of GHRT and 
diabetes incidence or O/E ratio was ob- 
served (P trend = 0.47 in the final Poisson 
regression model structure). In a similar 
manner, there was no statistical associa- 
tion between prescribed GH dose at 1 year 
and later occurrence of diabetes (P = 0.74 
in the final Poisson regression model 
structure). It is noteworthy that although 
diabetes incidence increased as expected 
with advancing attained age, the O/E ra- 
tio, although increased in the younger pa- 
tients, decreased with advancing attained 
age. Neither age at GHRT start, the num- 
ber of additional pituitary deficiencies, 
nor family history of diabetes had an im- 
pact on development of diabetes (data not 
shown). In addition, O/E diabetes ratio 
was not statistically significantly different 
in patients receiving glucocorticoid, 
males receiving testosterone replacement 
therapy, or patients who had undergone 
pituitary surgery prior to entry into KIMS 
compared with those who had no substi- 
tution therapy or surgery. 

In the subgroup of patients with 
available yearly fasting glucose and 
HbA lc data over a period of 6 years (239 
and 451 patients, respectively), plasma 
glucose concentration increased continu- 
ously from 84.4 ± 0.9 mg/dL to 89.5 ± 
0.8 mg/dL and HbAi c increased from 
4.74 ± 0.04% to 5.09 ± 0.13% (with 
the effect of age removed). This rep- 
resents a statistically significant increase 
of 0.70 mg/dL/year in plasma glucose 
concentration (95% CI 0.44-0.95; 
P trend < 0.0001) and of 0.036%/year 
in HbA lc value (0.024-0.048; P trend 
< 0.0001). 

CONCLUSIONS— This is the first 
study to report on the incidence of di- 
abetes mellitus in a large cohort of adult- 
onset GH-deficient patients receiving 
GHRT followed over a prolonged ob- 
servation period with a mean of almost 
4 years. The exclusion of patients with 
childhood-onset GH deficiency, Cushing 



care . diabe tesj oumals . org 



Diabetes Care, volume 35, January 2012 



59 



Diabetes incidence during GHRT 



Table 2 — Likelihood ratio statistics for type 3 analysis of the multiple Poisson regression 
model giving the magnitude of the explanatory value of each variable on the O/E 
standardized incidence ratio; reference population Kronoberg county, Sweden, 1998—2001 



Variahlp*t 

V (XL IttUlt 4> 


df 


2 

X 


Pr > v 2 


% v 2 
lo x 


% x 2 /df 


BMI (kg/m 2 ) 


i 


133.85 


<0.0001 


45.1 


45.1 


Year since GHTS 


i 


63.29 


<0.0001 


21.3 


21.3 


Year from diagnosis of pituitary 
















4.29 


0.0383 


1.5 


1.5 


Sex 


1 


5.00 


0.0253 


1.7 


1.7 


X 2 totals 




296.59 




100 





GHTS, GH treatment start; Pr, probability. *Additional deficiencies, family history of diabetes, calendar 
period, age at KIMS entry, age at diagnosis of pituitary disease, IGF-I SDS, and GH dose did not have any 
significant explanatory value in the multiple regression analysis. $A11 variables were numerical. (The difference 
in deviance between the fully categorized model and the numerical model was statistically nonsignificant.) 



disease, or acromegaly and those who had 
ever received GH ensured that confounding 
factors related to the cause or treatment of 
GH deficiency and possibly contributing 
to the development/manifestation of dia- 
betes were eliminated to the largest possi- 
ble extent. 

We have shown, using a sample of the 
Swedish general population as reference 
population (15), the O/E ratio of 6.0 
appeared significantly elevated in GH- 
deficient patients. This was also true 
when diabetes incidence in KIMS patients 
was compared with other reference popu- 
lations from different European regions 
and the U.S. (16-19); age- and sex-adjusted 
O/E ratios ranged from 2. 1 1 to 5.22. How- 
ever, the GH-deficient patients exhibited 
at baseline a substantial adverse risk pro- 
file for developing diabetes. Indeed, the 
patients developing diabetes were signifi- 
cantly older; had significantly higher BMI, 
waist circumference, triglyceride concen- 
trations, and blood pressure; and had 
significantly lower HDL-cholesterol con- 
centrations when compared with those 
who did not develop diabetes. It has pre- 
viously been shown that the presence of 



any of these risk factors considerably in- 
creases the risk for developing diabetes in 
the general population: 2-fold in the case 
of dyslipidemia or hypertension, 3-fold 
when BMI is 25-30 kg/m 2 , and 10-fold 
when BMI is >30 kg/m 2 (16). 

Because diabetes incidence varies in 
different regions of the world, O/E ratios 
were determined with reference popula- 
tions from not only northern Europe (15) 
but also central (17,19) and southern 
Europe (16) and the U.S. (18). And again, 
an increased age- and sex-adjusted diabe- 
tes incidence was documented in com- 
parison with these cohorts. All of these 
reports were population-based studies. 
The methods of ascertainment of dia- 
betes differed in the different studies: 
oral glucose tolerance tests (Bruneck and 
KORA S4/F4), fasting glucose concen- 
trations (Kronoberg), self-reported diagno- 
sis (U.S. adults and MONICA Augsburg), 
and any of the above-mentioned meth- 
ods plus an HbA lc value of ^6.5% or 
the use of any antidiabetes drug in the 
KIMS study. In addition to the higher 
prevalence of risk factors for diabetes in 
GH-deficient patients, these different 



procedures certainly influence the reported 
incidence rates. Self-reporting indeed will 
give lower incidence rates than studies re- 
lying on repeated testing (fasting blood glu- 
cose or oral glucose tolerance test) of a 
representative sample of the whole popu- 
lation of an area. In the Kronoberg study, 
opportunistic screening (i.e., testing all in- 
dividuals at routine or acute contact with 
the medical system, irrespective of the rea- 
son of consultation) was chosen. By includ- 
ing HbA lc and allowing all other possible 
ways of establishing the diagnosis of diabe- 
tes, we aimed at the most comprehensive 
way to identify manifestation of diabetes 
in any KIMS patient. We are well aware 
of the limitation of these different ap- 
proaches; however, the ideal age-, sex-, 
and risk factor-matched GH-deficient 
group not receiving GHRT in a prospec- 
tive controlled trial is not available. As a 
result of the established indication for 
GHRT in GH deficiency, it also is unlikely 
that such a study will ever be performed. 

Diabetes incidence was not associ- 
ated with the GH dose or with the IGF-I 
concentrations reached after 1 year, and 
diabetes incidence decreased with increas- 
ing duration of GHRT. These observations 
might be interpreted either as an indica- 
tion that GH does not play a causative role 
in the development of diabetes or, alter- 
natively, that diabetes develops above a 
GH threshold dose and IGF-I concentra- 
tions and that patients at risk for develop- 
ing diabetes do so early after initiation of 
GHRT. It is noteworthy that an increased 
incidence of diabetes also was observed in 
GH-deficient patients not receiving GH 
(2). Another observation deserves special 
attention. Although diabetes incidence in- 
creased with advancing (attained) age in 
GH-deficient patients as expected from 
studies in non-GH deficient populations, 
the O/E ratio decreased with advancing 
age. Additional hormonal deficiencies 
and their replacement therapies were sim- 
ilar in patients who developed diabetes 
and those who did not, as well as the prev- 
alence of a family history for diabetes and 
the age at diagnosis of pituitary disease. 
Also, O/E diabetes ratios were similar in 
patients receiving glucocorticoid, males 
receiving testosterone replacement ther- 
apy, and patients who had undergone pi- 
tuitary surgery prior to KIMS entry when 
compared with those not receiving these 
therapies or without prior surgery. 

The mechanism of the increased di- 
abetes incidence in GH-deficient patients, 
whether or not on GHRT, is not fully 
understood. In GH-deficient patients not 



Table 3 — Contrast estimates from the final multiple Poisson regression model 



Effect on* 


Estimate* 


95% CI 


O/E per year from pituitary disease to KIMS entry 


1.0136 


(1.0009-1.0266) 


O/E per year since first GH treatment (since KIMS entry) 


0.8467 


(0.8102-0.8848) 


O/E per unit BMI (kg/m 2 ) 


1.0938 


(1.0783-1.1095) 


O/E sex ratio: female vs. male 


1.2290 


(1.0261-1.4720) 



*0/E decreased by —4% per year of attained age (P < 0.0001), with somewhat lower decrease at attained 
ages >50 and somewhat greater decrease at lower attained ages {background attended age model: Model exp 
[-0.0393 X (age - 50) + 0.001 X (age - 50) 2 1). The model estimates that O/E 1) increases by 1.36% peryear 
for year from diagnosis of pituitary disease to KIMS entry (95% CI 0-2.66), 2) decreases by 15.3% by year for 
year since GH treatment start (11.5-19.0), 3) increases by 9.4% per unit BMI (7.8-10.95), and 4) is 22.9% 
higher in females than in males (2.6-47.2). 
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receiving GHRT, the adverse risk profile 
for diabetes can be considered an impor- 
tant contributing factor (21). This adverse 
risk profile is positively influenced by 
GHRT (21), but GH also opposes insulin 
action. Several studies in non-GH defi- 
cient control subjects and GH-deficient 
patients suggest a complex role of GH 
in carbohydrate metabolism, with free 
fatty acids playing an important role 
(1,22,23). 

Another important finding of the 
current study is the small but significant 
increase of fasting glucose concentrations 
of 0.7 mg/dL/year and HbA lc levels of 
0.036%/year. Although this seems low 
at first glance, it represents an increase 
of fasting glucose of 7.0 mg/dL and of 
HbA lc of 0.36% in 10 years, changes 
that must be considered as clinically rel- 
evant. However, these figures need to 
be compared with observations in GH- 
deficient patients without GHRT and in 
control subjects. In a recent smaller study, 
comparable increases in fasting blood glu- 
cose and HbA lc concentrations were ob- 
served in GH-deficient patients receiving 
GHRT and those not receiving GHRT 
(24). The increase in blood glucose in 
that study, with a mean observation pe- 
riod of 31.7 months, was 7.4 mg/dL (2.84 
mg/dL/year) in untreated patients, a value 
considerably higher than the one observed 
in our cohort. 

In summary, incidence of diabetes 
irrespective of reference population ap- 
pears to be increased in this study of 
5,143 patients followed over a mean 
observation period of 3.9 years. Patients 
who developed diabetes were signifi- 
cantly older, had a more adverse lipid 
profile, and had higher BMI and blood 
pressure at baseline. Fasting blood glucose 
concentrations and HbA lc also increased 
significantly; however, in another recent 
study comparing GH-treated to untreated 
GH-deficient subjects, no difference was 
observed. Thus, despite widely demon- 
strated benefits of GHRT in adult GH de- 
ficiency (25), a marked increase in the risk 
of developing diabetes must be consid- 
ered. All patients, but particularly those 
with an adverse risk profile, should be 
carefully followed regarding parame- 
ters of glucose metabolism. 
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